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ABSTRACT 


Multispectral thermal 1R data were acquired with NASA's Bendlx 24-channel 
scanner flown in an aircraft over the E. Tlntic Utah mining dist <ct. These 
digital image data required extensive computer processing in order to put the 
information into a format useful for a geologic photo-interpreter. Simple 
enhaicemer.t procedures were not sufficient to reveal the total information 
content because the data were highly correlated in all channels. The data 
were shown to be dominated by temperature Variations across the scene, while 
we were interested in the much more subtle spectral variations between the 
different rock types. In this paper we discuss the image processing tech- 
niques employed to analyze these data. 

The first step was to create a view-angle-corrected image from each tape 
to evaluate data quality and select the area of interest. Next, following the 
work of Vincent and Thomson (1972), ratio images were produced for each pair 
of adjacent wavelength channels, and a color ratio composite created. To 
quantify the data further, on-board sensor calibration and an atmospheric 
model utilizing Salt Lake City radiosonde data were used to determine the 
surface brightness temperature from the 11-12 pm channel. This surface 
brightness temperature, combined with the calibration and atmospheric model, 
was then used to create emissivity images for each of the remaining wavelength 
channels. A color ratio composite of three emissivity images provided im- 
proved separation but uas not fully satisfactory, due in part to increased 
noise. Finally a principal component transformation, followed by a Gaussian 
stretch, followed by an inverse transformation to the original axes achieved 
satisfactory separation of many of the major rock types. This procedure was 
invoked for both the raw radiance data and the emissivity images. While both 

v 
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revealed about the same amount of rock- type information, the former was Judged 
to be much easier to interpret because of the topographic information intrin- 
sic in the temperature was still contained in the final product. 


INTRODUCTION 


Laboratory measurements of middle-infrared (MI r , 5-MO pm) spectra of 
rocks and minerals show numerous diagnostic features. The region between 8 
and 14 pm holds the most promise for remote sensing because this is an excel- 
lent atmospheric window and also the region of maximum thermal emission at 
terrestrial surface temperatures. Within this spectral range, the most prom- 
inant spectral features a^ise from the silicon-oxygen stretching vibrations. 
These features change location and intensity with varying composition and 
structure (Hunt and Salisbury, 1974, 1975 t 1976; Vincent et al . . 1975). Lyon 
(1965) indicated that these intense features (known as reststrahlen bands in 
reflectance spectra) shift to shorter wavelengths with increasing silica 
content. 

The possibility of exploiting these spectral features for remote sensing 
of rock-type from aircraft or satellite has been suggested by numerous authors 
(Vickers and Lyon, 1967; Vincent and Thomson, 1972; Vincent, 1975). However, 
owing to lack of multispectral scanners (Vincent, 1975), very few tests of the 
technique have been possible. Hovis et al . (1968) flew a non-imaging emit- 
tance spectrometer over desert terrain in California. They concluded that 
even though atmospheric effects were significant, the reststrahlen bands of 
silicates were o&servable. Lyon (1972) flew a non-imaging spectrometer over 
the Pisgah lava flow3 and surroundings. He was able to distinguish between 
the basalts, the alluvial fans, and the windblown sands, and also recognized 
spectral differences within the basalt flows. Two tests with a two-channel 



images on which they could distinguish between the quartz sand or sandstone, 
and the non-silicate surface material. Vincent and Thomson (1972) flying a 
scanner with the same two channels over the Pisgah Crater area were able to 
discriminate dacite from uasalt, and a rhyolitic tuff was distinguished from 
the surrounding alluvium. However, until recently, lack of adequate imaging 
multispectral scanners has precluded further validation of the use of spectral 
emittance data. 

In November of 1975, NASA’s 24-channel scanner was flown over the East 
Untie mining district of central Utah. On this flight six channels of radi- 
ance data were acquired in the 8.27 to 13 pm range. This flight provided an 
unprecedented opportunity to analyze multispectral mid- T R image data in an 
area of moderate geologic complexity. The district, in an area of high relief 
and moderate vegetation, consists mainly of Tertiary silicic igneous rocks and 
Paleozoic quartzite and carbonate rocks which have been locally hydrothermal ly 
altered. 


Improved image processing algorithms were required to best utilize these 
data sets. Simple contrast enhancement and color compositing of the various 
spectral channels wer^ insufficient to reveal the information content, because 
the data were highly correlated in all channels. The data were dominated by 
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temperature variations across the scene, while we were interested in the much 
more subtle spectral variation:: between the different rock types. This paper 
describes the image processing applied to these data sets, with illustrations 



2 


IMAGE PROCESSING 


I 


\ 





✓ 


A. Geometric Corrections 

The raw data used in this study were obtained by NASA's Bendix 24-channel 
scanner. Three parallel flight lines were flown, with the aircraft at an 
elevation of between 5000 m and 6000 m above the terrain. The total of 12 
computer compatible tapes (CCT's) were produced at Johnson Space Center, with 
each tape containing all 24 channels (operative or not) for a portion of a 
flightline. Each picture element on the original data covered a ground area 
of approximately 10 by 10 m at the sub-aircraft point. 

Because the recorded data values represented radiance measurements for 
equal angular displacements of the scanning apparatus, a panorama correction 
was applied immediately after logging the data. To accomplish this correc- 
tion, the raw data were resampled along each scan line such that each data 
value represented the same surface area, as described by Gillespie and Kahle 


(1977). 

B. Contrast Enhancement and Color Display 

Of the 24 channels, there were six operational channels in the mid-IR 
wavelength. Table I shows the position and effective bandwidth3 in the 
electromagnetic spectrum of these channels. 


Table I 

Mid-IR Channels on 24-Channel Scanner 
Channel Bandlimits (micrometers) 


17 

18 

19 

20 
21 
22 


8.3 - 8.8 
8.8 - 9.3 

9.3 - 9.8 
10.1 - 11.0 
11.0 - 12.0 
12.0 - 13.0 


3 








The images of the ground scene for these six channels were Individually 
contrast enhanced, and those from one tape near the center of the test area 
are displayed in Figure 1. The contrast enhancement used here was obtained by 
matching the cumulative distribution functions of each of the six channels to 
a reference truncated Gaussian distribution of + 2.0 <T. 

It is immediately apparent in Figure 1 that surface temperature, control- 
led by topography, is the dominant information in these images. Only very 
minor changes can be detected between images. Channel 22 has noticeable noise 
striping. 

In order to better perceive subtle brightness changes in the scene, it is 
customary to combine three of the images into a single color display. This is 
accomplished by assigning each image a primary color and simultaneously pro- 
ject iiig all three images onto color film. Figure 2 illustrates this process 
using channels 17, 18, and 20 as blue, green, and red respectively. The 
original radiance data had some areas that were saturated, i.e., the ground 
temperature was either too hot or too cold for the measured radiance to lie 
within the range of the recording system. While these dices apj eared eitner 
white or black in the data displayed in Figure 1, we have set all these areas 
to black in the color composite, Figure 2, and have done so in most subsequent 
images. 

In order to compare directly the processed images with geologic informa- 
tion, a registration procedure was carried out using the Universal Transverse 
Mercator (UTM) grid system. First the image displayed in Figure 1 was conca- 
tenated with the image of the area lying Just to the north. Tiepoints were 
then chosen that could be identified both on this large image and on the USGS 
topographic map of the area. On the image, the tiepoints were referenced by 
their line and sample coordinates and on the map by their UTM coordinates. 


4 









W-' t 





Tho set of tlepoints was inert fcc determine a second order polynomial transfor- 
mation which mapped the line and sample tiepoint coordinates to their respec- 
tive UTM coordinates. Using this polyncminal transformation, » regular tie- 
point grid was generated that covered the entire image. The displacements 
associated with the grid points then controlled a geometric transformation 
algorithm by which the full image was resampled. Figure 2 and all subsequent 
color images display the results of this geometric transformation. 

In Figure 2 we note that the dominance of temperature effects, with the 
resulting high inter-Dand correlation, causes the range of color to be quite 
low. Nevertheless it is possible to recognize the location of the Tintic 
Quartzite — the distinctive pink unit along the left edge of the image. 

Beyond that distinction we see only that the image is divided into light pink 
or gray green. 

C. Band Ratioing 

One technique used to deal with the problem of high correlation between 
channels is to create ratios of one channel to another ( Vincent and Thomson. 
1972; Vincent, 1973). In this way, the temperature effects tend to cancel 
out. However, became of the high inter -oand correlation, the resultant ratio 
images have very narrow distributions of data values centered about unity. 

When these ratio images are contrast enhanced for display purposes, high gains 
must be applied in order to produce images utilizing the full black to white 
brightness range. Because random noise is intrinsically uncorrelated from 
channel to channel , the above procedure reduces the apparent signal-to-noiae 
ratio. Consequently, the application of large gains to the ratio images tends 
to emphasize the uncorrelated no.'se, as car. be seen in Figure 3 which displays 
the ratio images of adj? lent bands: 17/18, 18/19, 19/20, 20/21, and 21/22. 

As before, it is useful to construct a single color display, this time using 
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three ratio images. Figure 4 displays the color ratio composite using 17/18, 
18/19, 19/20, as blue, green and red respectively. In this image more color 
separation has become apparent. The Tintic Quartzite along with other quartz- 
rich rocks are a deep purple. However, the rust colored areas and green areas 
each consist of several rock types, alluvial material, and/o:' vegetation. The 
magenta appears to be temperature related, appearing primarily on cool north- 
facing slopes. 

D. Sensor Calibration, Atmospheric Corrections and Derivation of Emittance 

The black and white images above, together with the associated color 
composites represent the initial products that were created. The existence of 
onboard calibration equipment made it possible to generate a -physically mean- 
ingful quantity, the surface spectral emittance, which potentially might be 
more directly interpretable. We use the term emittance, rather than emissiv- 
ity, to indicate that we are considering a surface in its natural state, 
rather than an ideal specimen. To determine the emittance value corresponding 
to a given raw data value required a simple atmospheric model which provided 
average values in each channel of atmospheric transmission, emission, and 
scattering properties. Using the onboard calibration and the atmospheric 
model, it is possible, as described below, to determine the surface tempera- 
ture at every data point if we assume that the emittance in the longest wave- 
length channel (12 - 13 ^m) is constant across the entire image (Lyon, 1965; 
Vincent, 1972). Once this temperature is known, it is then possible to calcu- 
late the surface emissivity for every point in the scene for each of the 
remaining wavelength channels. 

Following the notation of Vincent and Thomson (p. 2467) the total spec- 
tral radiance at the sensor for the wavelength \ is given by 
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l8 xHbX (T) + (1 - e x> L skyx ] T AA + L VX 


( 1 ) 



L 


where 


is the surface emittance, 

is the spectral radiance of a black body at surface 
temperature T, 

L , is the spectral radiance incident upon the surface, from the 

sV.y>. 

atmosphere, 

t, is the spectral atmospheric transmissivity, 

AX 

L is the spectral radiance from atmospheric emission and 

V A 

scattering that reaches the sensor. 

In this study, equation (1) is utilized in an average sense. That is, it is 
written once for each channel , using the average values from the atmosphere 
model for L 9k yX» and L y ^. The intent is to be able to solve the equation 

for the average emittance at each data point in the image for all the thermal 
channels. For a given data point, these average emittances can then be 
plotted as a function of wavelength. The atmospheric model employed was 
LOWTRAN 4 (Selby, et al .. 1978). Salt Lake City radiosonde data for the day 
of the overflight was used in the model to define atmospheric conditions. The 
model then produced atmospheric transmissivity, x^, atmospheric directional 
upward radiance, L , and directional downward radiance. The latter was 

^ A 

integrated over the sky hemisphere to determine the downward radiance at the 



surface, L^. 


These quantities are given in Table II. 
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Table II 


Atmospheric Parameters from LOWTRAN 4 


Wavelength 


Interval 

t AA 

L skyA 

L vA 

( pm) 


2 

(watts/m sr urn) 

2 

(watts/m sr pm) 

8.3 - 8.8 

0.848 

4.198 

0.962 

8.8 - 9.3 

0.885 

3.286 

0.650 

9.3 - 9.8 

0.801 

5.116 

1.180 

10.1 - 11.0 

0.910 

3.248 

0.603 

11.0 - 12.0 

0.912 

2.696 

0.498 

12.0 - 13-0 

0.816 

4.945 

1.143 


Referring back to equation (1), once the atmospheric parameters are 
defined, the remaining unknowns are the surface emittance, and the surface 
temperature T. L^.^T) can be calculated using Planck's blackbody radiation 
function once T is known. The following derivation illustrates how a corres- 
ponding surface temperature linage may be constructed using the calibration 
information provided and one of the thermal channels. 

For a given thermal channel i, the relative sensor response R*(A) has a 
Gaussian-like form with response falling to zero at the extremes A* and A* 
(i.e. R i (A) = 0 for A < a| and A > A*). The wavelength intervals given in 
Table II are the approximate domains for which R i (x) is greater than half its 
maximum value. 


To deal with the non-uniform sensor response, all spectral radiance 
functions in equation (1) are weighted by the relative filter response R*(A) 
and integrated over the bandpass to obtain an effective radiance value 




' 7 


R'vA) L (A) dA 


( 2 ) 


where the overbar will be used '.o denote effective radiance. The normaliza- 


tion constant N 1 is chosen to make the effective radiance equal to that value 
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of constant spectral radiance which would produce an equivalent respo. je in 
the sensor. This constant is then 


N 


i 



Aa) dA 


The response of the sensor in teros of output data value (DV) is then assumed 
to be a linear function of this effective radiance according to 


DV 


i 


A. L 1 + B. 
i i 


(3) 


The radiance of a blackbody source at temperature T is given by 


■4 ( X - T) 


ttA" 


c./AT 

(e Z - 1) 


watt 

2 

m sr 


where 

c.j = first radiation constant = 3*7'4 x 10~ ° watt m , 
c 2 = second radiation constant = 0.0144 m°K. 

Then the effective or normalized radiance from a blackbody source in the same 
wavelength region is given by applying equation (2) 


(T) 




c 1 R i (A)dA 

c./Xt ' 
(e 


1) 


(4) 


The calibration system consisted of two blackbody sources maintained at 
different, known temperatures. After the completion of each scan line of the 
ground scene, the two blackbodies were scanned and the corresponding data 
values were recorded along with temperature. Let the temperature of the 
colder blackbody be T^ and that of the hotter blackbody be T H> The two cali- 
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bration points then became 


K- ^b «c»- K> ^b <V> 


where the blackbody effective radiances E bb (T c ) and L bb (T H ) are calculated 
according to equation (4). Using the linear sensor response model of equation 
(3) i a system of two equations is generated 


DV J ’ A i <V + B i 

DV J ' A i ^b l V + B i 


which can easily be solved for the coefficients A.^ and 

When the scanner samples the actual surface, all the radiance terms of 
equation (1) influence sensor response, so the linear sensor model becomes 


i 


DV 


A i l L 


(T.e^ + B. 


(5) 


where 




\ ^b (T) + a ' C l ) \ B sky + 


( 6 ) 


The latter expression is obtained by applying equation (2) to the right side 
of equation (1). The blackbody effective radiance L bt) (T) is determined 
according to equation (4). The other two effective radiance terms are easily 
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relatel to the values in Table II by assuming L 


s«cyA 40(1 KrX £re con3tant 


through the wavelength region. Thus 


* 4 , 


n 

i I R 1 (A) L 1 . . 

N 1 j , Sky> 


dA 


and similarly, 



From laboratory spectra, Vincent et al . , 1975 state that in the 12 - 13 pm 
region, the emittance of many silicate igneous rocks approaches a constant 
values £*0.96). This wavelength region corresponds to channel 22. If this 
result is applied to equation (5): 

Dv22 ’ A 22 E “t <T ' £ 22 ' °- 96 > + B 22 (7> 

We now have an expression involving only data values and temperature. Since 
the temperature T cannot be solved for analytically, a numerical inversion 
technique was employed. The signal level, DV“^, was calculated from equation 
(7) for a sufficiently broad range of temperatures to encompass the allowable 
range of DV. This information was stored in tabular form. The table was then 
inverted by interpolation, yielding a tenperature value for each possible 
value of DV^. The table was then referenced for each picture element to 
produce a temperature image. 
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The emittance values for channels 17 - 21 could then be determined in the 


I 


h 

s 


following manner. Suppose that the constants and were determined as before. 
Solving equation (6) for L 1 (T,e^): 


l Lt <T - e i> 


i- (DV 1 - B,) 


Substituting into equation (A) and solving for gives: 


^ (T) * i - E ; 


ri 

sky 


T- (DV 1 - B.) 
A i 1 


L 1 - T 1 * L 1 

v A sky 




I 



Now, at each data point in the image, T has already been determined and DV* is 
known. So, the emittance at that data point is completely determined. In 
this manner, the corresponding emittance images for each of channels 17 
through 21 can be determined using the temperature image created from channel 
22 along with the calibration information. (Note again that the emittance for 
channel 22 is assumed to be constant.) 

With the available data, channel 22 was of such poor quality that using 
it to generate the temperature image resulted in unacceptably noisy emittance 
images. Channel 21 was used instead, with €21 * 0.93, but assuming a constant 
emittance for this channel is questionable. In this wavelength region (11 - 
12 um) many of the silicates have a steeply rising spectrum, and the carbon- 
ates have a narrow but intense band (Hunt and Salisbury, 1974, 1975, 1976). 
Nevertheless, use of channel 21 was considered to be preferable to use of 
channel 22 with its excessive noise. The temperature image based on this 
calculation is shown in Figure 5. Figure 6 displays contrast enhanced emit- 
tance images for channels 17, 18, 19, and 20, derived using channel 21 to 
generate the temperature image. Figure 7 is a color composite formed using 


A. 
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a. Channel 17 


Channel 18 


c. Channel 19 


Figure 6. Emittance images for channel 17, 18, 19, and 20 using temperature 
from Figure 5. 







channels 17, 18. and 20 of the above emlttance Images displayed as blue, 
green, and red respectively. 

In the Individual channel emlttance images, it is apparent that the 
quartzite (left center) has low relative emlttance in all wavelength channels, 
although not as pronounced in the highest wavelength. Discrimination of the 
rock units does not appear to be based upon the moving of the reststrahlen 
band to higher wavelength with decreasing silica content so much as just the 
existence of a very broad intense spectral feature across all wavelengths in 
the high quartz content rocks. In the color composite emlttance image quite a 
bit of rock type discrimination has become possible, as will be discussed 
further with Figures 8 and '. 

Several other atmospheric models and several values of con3tant emiss»iv- 
ity for channel 21 were used in the foregoing derivation of emissivity 
images. The changes in constant emissivity value had no noticeable effect on 
the emissivity images — merely raising or Covering the absolute value of the 
derived emissivity values. Evaluation of the emissivity images derived using 
different atmospheres show minor but not significant changes. Ability to 
distinguish geologic units was unchanged. 

E. Color Enhancement 

As is apparent, particularly in Figure 2, but also in Figures 4 and 7, a 
problem associated with highly correlated multispectral data is a lack of 
color range in a color composite display. This problem was solved for the 
emlttance and thermal radiance composites by use of a multispectral transfor- 
mation technique that removes correlation between components without rotating 
the data distribution. 

The simplest and most frequently used approach to enhancement of a color 
image is to separately perform contrast enhancement on each individual image 
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Figure 8. Color enhanced emittance composite. 
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component. To produce Figure 2, for example, individual intensity transfor- 


mations 


T i (x i } 


were performed for the i th spectral components (i = 17, 18, 20), where x^ 
represents the value of an arbitrary picture element in the i th component 
before enhancement, and Y^ its value after enhancement. For Figure 2, Y^ is 
displayed as blue, Y^ as green and Y^q as red. For scenes with high inter- 
band correlation, seDarate intensity transformations merely lengthen the 
narrow extended histogram distribution of the original, without accomplishing 
the high degree of spreading required to use the full available range of color 
to display differences in the scene. This situation is illustrated by the 
two-dimensional histogram of Figure 10. To achieve the needed spreading, joint 
transformations of the form 

(Y,y» Yiq , Y2q) = D (x-|y, *ig» x 2Q^ (®) 


are required. 

The technique used to accomplish this joint transformation in the present 
application involved several steps (Soha and Schwartz, 1978). First, a prin- 
cipal component transformation was applied to the three components (17, 18 and 
20) of the emittance image. If we regard the original three components as 
forming a cartesian 3-space, then the principal component transformation 
consists of a rotation of the coordinate axes, producing new components which 
are linear combinations of the original. The transformation is implemented 
via a 3 x 3 rotation matrix applied to each multispectral picture element 
vector (x-jy. x.^, '^ q j to P rod uce three new component images, according to 
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(z i , z 2 , ) 


P (x^* x l8 , x 2Q ) 


The matrix is derived from the 3x3 covariance matrix C, which is estimated 


over a statistically significant subset of the image by 


Z (x ± - x ± > (x. - x.) 


where i, j = 17, 18, 20 and the summation is over N picture elements and x^ is 
the mean value of component i. The diagonal elements of C (with i = J) are 
the variances of the three components; the off-diagonal elements (i i j) 
describe the correlation between components. For highly correlated data, as 
in Figures 2 4 and 7, these cff-diagonal elements have large absolute 

values. T he difficulty in displaying highly correlated data can be thought of 
as inefficient use of available display dynamic range: "duplicate" informa- 

tion is displayed in more than one channel, with a rasulting waste of poten- 
tial discriminating power. Fortunately, it is possible to rotate the coordin- 
ates to achieve new uncorrelated components. After this rotation P, the 
covariance matrix Cp relating the new components has no nonzero off-diagonal 
elements. This new diagonal covariance matrix Cp derives from the original 
covariance matrix C by a similarity transformation 


PCP 1 = C c 


Since Cp is diagonal, both it and the rotation matrix P can be found by solv- 
ing the eigenvector equation 


33 


PC 


(Further description of uses of the principal component transformation can be 
found in Ready and Wintz, 1973.) Production of new uncorrelated components 
using the principal component transformation P was the first step in the color 
enhancement procedure. 

Next, each rotated component was contrast enhanced, using standard indi- 
vidual intensity transformations. Because the new components are uncorrela- 
ted, individual contrast transformations were now capable of producing the 
desired spreading in the histogram. All components were enhanced (using 
individual non-linear intensity transformations, collectively referred to as G 
below) so that their histograms approximated Gaussian distributions of a 
specified variance, resulting in a three-dimensional composite histogram which 
approximates a spherically symmetric three-dimensional Gaussian and hence 
occupies the majority of the available color space. 

After the variances of the principal component images have been equa- 
lized, it is possible to rotate the "multispectral" coordinates further with- 
out re-introducing correlation (Taylor, 1973 ). To produce the final color 
enhanced image, a coordinate rotation which is the inverse of the original 
principal component rotation was applied. Thus the entire joint transforma- 
tion of equation (8) is defined by 

D z P" 1 * G * P 

The resulting color enhanced product is shown in Figure 8. The inverse rota- 
tion can be viewed as returning to the original coordinate orientation. This 
operation has the advantage of restoring basic color relationships which 
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existed in the original scene and hence facilitating certain Inferences. For 
example, red areas in Figure 8 result from proportionately lower emittance 
values in the channel 17 wavelength band. 

The same color enhancement procedure was applied to the radiance compo- 
site (the data in Figure 2), with the resulting product shown in Figure 9. 

This product has proven to be preferable for interpretation, primarily because 
it provides visual separation of topographic and emissivicy variations. The 
first principal component image is a weighted average of the three input 
thermal radiance channels, and is in fact almost entirely dependent upon 
ground temperature. When the correlation between the first principal compo- 
nent image and the temperature image of Figure 5 is computed, a correlation 
coefficient of .967 is found. On the other hand, when the first principal 
component image is compared to the emissivity images for the channel 17, 18 
and 20 wavelength bands, correlation coefficients of less than 0.03 result, 
indicating little dependence of the first principal component upon emissiv- 
ity. After the inverse rotation step of the color enhancement procedure, the 
first principal component, which was very nearly the average of the original 
components, becomes essentially tne intensity component of the final color 
composite. Since the ground temperature variation of Figure 5 is largely 
topographically controlled, intensity variation in the final product of Figure 
9 results primarily from topography. When the second and third principal 
component images are compared with the temperature image, minimal correlation 
coefficients (.038 and .032) result. These lower order principal components 
in fact depend primarily upon emissivity differences. After the inverse 
rotation, variations in these components become color variations in the final 
product. Thus in Figure 9 color differences are indl" !5 Hve of emissivity 
differences in ground materials in the scene. These colors are modulated by 
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intensity variations which reflect topography, and hence are useful for estab- 
lishing geographic orientation during interpretation. By comparison, in 
Figure 8 both color and intensity variations indicate emisslvity differences. 

In Figure 8, but particularly in Figure 9, we are now able to distin- 
guish several rock types. Field checking of the image has revealed that most 
of the distinct color boundaries correspond to rock-type boundaries, although 
not all rock-type boundaries are distinguishable. In Figure 9, the color red 
is a good indicator of quartz content — the more intense and "true 1 * the red, 
the higher the quartz content. In general, the pink areas are silicified 
areas, sandy limestones, quartz latites, and quartz monzonites; the blues 
correspond to latites and monzonites and some carbonates ; the purples to 
clays; and the green areas are carbonates, vegetation and alluvium. 

It is not the intent of this paper to describe the geological interpreta- 
tion of the images or their spectral significance. These are being described 
elsewhere. For instance, for a comparison of Figure 9 with a simplified 
geologic map, see Kahle and Rowan (1979). However, it is clear that with the 
image processing techniques used here, it is possible to separate several rock 
and ground-cover types based on their spectral emittance differences. 

It is unfortunate, with the great potential shown here for geologic 
remote sensing, that the 24-channel scanner is no longer operational nor does 
there exist any other scanner with more than two mid-IR channels. This work 
reiterates and emphasizes the need for development of a multispectral, mid-IR 
scanner system for aircraft and satellite remote sensing. 
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